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' We have analysed the data collected by OPAL at centre-of-mass energies between 189 and 
209 GeV searching for Higgs boson candidates from the process e"'"e~ — >■ h^Z'' followed by the 
^ . decay of h° —* A°A° where A*^ is the CP-odd Higgs boson. The search is done in the region 
where the A° mass, ttia, is below the production threshold for bb, and the CP-even Higgs 
boson mass m-h is within the range 45-86 GeV/c^. In this kinematic range, the decay of 
h° —>■ A°A° may be dominant and previous Higgs boson searches have very small sensitivities. 
This search can be interpreted within any model that predicts the existence of at least one 
scalar and one pseudoscalar Higgs boson. No excess of events is observed above the expected 
Standard Model backgrounds. Model-independent limits on the cross-section for the process 
e"'"e~ — i> h'^Z° are derived assuming 100% decays of the h° into A^A*^ and 100% decays of 
the A^A" into each of the following final states: cccc, gggg, t^t^t^t^ , ccgg, ggr"'"r~and 
ccr+r^. The results are also interpreted in the CP-conserving no-mixing MSSM scenario, 
where the region 45 < < 85 GeV/c^ and 2 < mA < 9.5 GeV/c^ is excluded. 



To be submitted to European Physics Journal C 



G.Abbiendi^ C.Ainsley^ P.F.Ak esson^, G. Alexander^^, J.Allison^^, P. Amaral^, 
G.Anagnostou^ K.J. Anderson^, S.Arcelli^, S.Asai^^, D.Axcn^'^, G. Azuelos^^'", I. Bailey^*^, 
E.Barberio^ R.J. Barlowl^ R.J. Batley^ RBechtle^^, T.Behnke^s, K.W.BelF, RJ.Bell^, 
G.Bella22, A. Bellerive^ G. Benelli^ S.Bethke^^, Q.BiebeF, I.J. Bloodworth^, O.Boeriu^°, 
P.Bock^^ D. Bonacorsi^, M. Boutemeur^^ S.Braibant^, L. Brigliadori^, R.M. Brown^°, 
K.Buesser^^, H.J. Burckhart^, S. Campana^, R.K. Carnegie^, B.Caron^^, A.A. Carter^^, 
J.R.Carter^ C.Y. Chang^^ D.G. Charltonl•^ A. Csilling^'S, M. Cuffiani^, S.Dado^^, 

G. M. Dallavalle^, S. Dallisoni^, A. De Roeck^ E.A. De Wolf«, K.Dcsch^^, B.Dienes^^^, 
M.Donkers^ J. Dubbert3^ E.DuchovnP, G. Duckeck3^ LP. Duerdoth^^ E. Elfgren^^ 
E.Etzion22, F.Fabbri^, L.Feld^°, P. Ferrari^ F. Fiedler=^\ I. Fleck^^, M. Ford^ A. Frey^ 
A.Furtjes^ P. Gagnon^^ J.W. Gary^ G.Gaycken^^, C. Geich-Gimbel^, G. Giacomelli^, 

P. Giacomelli^, M. Giunta^, J. Goldbcrg^^, E. Gross^"', J. Gmnhaus^^, M. Gruwe^, 
P.O.Giinther^, A. Gupta^ CHajdu^^, M. Hamann^s, G.G.Hanson^ K.Harder^s, 
A. HareP^ M. Harin-Dirac^ M. Hauschild^ J. Hauschildt^^, CM. Hawkes^ R. Hawkings^ 

R.J.HemingwayS, C.HenseP, G. Herten^o, R.D.Heuer^^, J.C.HilP, K. Hoffman^, 
R.J. Horner^, D. Horvath^^'"^, R. Howard^^, P. Hiintemeyer^^, P. Igo-Kemenes^^, K.Ishii^^, 

H. Jeremie^^, P. Jovanovic^, T.R. Junk^, N. Kanaya^^, J. Kanzaki^^, G. Karapetian^^, 

D.Karlen^, V. Kartvelishvili^^, K. Kawagoe^'\ T. Kawamoto^^, R.K. Keeler^^, 
R.G.Kellogg^^ B.W. Kennedy^^, D.H. Kim^^ K. KW^ A.Klier^^, S.Kluth^^, 
T. Kobayashi^^, M. Kobel^, S. Komamiya^^, L. Kormos^^, R.V. KowalewskP, T. Kramer^^, 
T. Kress^, P. Krieger^'^, J. von Krogh^^, D.Krop^^, K.Kruger^, M. Kupper^^, 
G.D. Lafferty^^, H. Landsman^\ D. Lanske^^, J.G. Laytcr'^, A. Leins^\ D. Lellouch^^, 
J.Lctts^^ L.Levinson24, J.Lillichi°, S.L. Lloydl^ F.K. Loebingerl^ J.Lu^^ J.Ludwig^", 
A. Macpherson^S'*, W.Mader^, S. Marcellim^, T.E. Marchant^^ A.J. Martin^^ J.P. Martin^^ 
G.Masetti^ T.Mashimo2^ P.Mattig™, W.J. McDonald2^ J.McKenna^^ T.J.McMahon^ 
R.A.McPherson^^, F. Meijers^, P. Mendez-Lorenzo^^ W. Menges^^, F.S. Merritt^, H.Mes^'", 
A.Michclini^, S. Mihara^^, G. Mikenberg^^, D.J.Millerl^ S.Moed^^, W.Mohr^o, T.MorP, 

A.Muttcri", K.Nagai^^^ I. Nakamura^^, H.A.NeaP, R.Nisius^^, S.W. 0'Neale\ A. Oh^, 
A.Okpara", M. J. Oreglia^ S. Orito^^, C.PahF, G. Pasztor^'^, J.R.Pater^^, G.N. Patrrfo, 
J.E.Pilcher^, J.Pinfold^^, D.E.Plane^ B.Poli^, J. Polok^ O. Pooth^^ M. Przybycien^'", 
A.Quadt^ K.Rabbertz^ C.Rembser^ P.RenkeP, H.Rick^, J.M.Roney^^, S. Rosati^ 
Y.Rozen^i, K.Runge^°, K. Sachs^, T. SaekP, O. Sahr3\ E.K.G. Sarkisyan^-^ 
A.D. Schaile^^, O.Schaile^^, P. Scharff-Hansen^, J. Schieck"^^, T. Schorner-Sadenius^, 
M.Schr5derS, M. Schumacher^, CSchwick^, W.G.Scott^o, R. Seuster^^'^, T.G. ShearsS''^, 
B.C.Shen^ C.H. Shepherd-Themistocleous^ P. Sherwood^^ G.Siroh^ A.Skuja^^ 

A. M.Smith^ R.Sobie^e, S. Soldner-Remboldi°'<^, S. Spagnolo^o, F.Spano^ A.Stahl^, 
K. Stephens^^, D.Strom^^, R. Strohmcr^^ , S.Tarcm^^, M. Tascvsky^, R.J. Taylor^^, 

R. Teuschcr'^ , M. A. Thomson^, E. Torrence^^, D.Toya^'^, P. Tran"^, T. Trcfzger^^, A.Tricoh^, 

I. Trigger^ Z. Tr6csanyP'% E. Tsur^^, M.F. Turner- Watson^ I. Ueda2^ B. UjvarP■^ 

B. Vachon^*^, C.F. Vollmer^\ P. Vannerem^°, M. Verzocchi^^, H.Voss^ J. Vossebeld^''^, 
D.Waller^, C.P.Ward^ D.R.Ward^ P.M.Watkins^ A.T.Watson^, N.K.Watson^, 

P.S.Wells^ T.Wengler^ N.Wcrmes^, D. WetterUng^^ G.W. Wilsonl6'^ J.A. Wilson^ 
G.WolP, T.R.Wyattl^ S. Yamashita^^, D.Zer-Zion^ L.Zivkovic^^ 



^School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
^Dipartimento di Fisica dell' Universita di Bologna and INFN, 1-40126 Bologna, Italy 



1 



^Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
^Department of Physics, University of California, Riverside CA 92521, USA 
^Cavendish Laboratory, Cambridge CBS OHE, UK 

^Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 
Ontario KIS 5B6, Canada 

^CERN, European Organisation for Nuclear Research, CH-1211 Geneva 23, Switzerland 
^Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 

USA 

^°Fakultat fiir Physik, Alb ert-Ludwigs- Universitat Freiburg, D-79104 Freiburg, Germany 
^^Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 
^^Indiana University, Department of Physics, Bloomington IN 47405, USA 
^^Queen Mary and Westfield College, University of London, London El 4NS, UK 
^^Technische Hochschule Aachen, HI Physikalisches Institut, Sommerfeldstrasse 26-28, D- 
52056 Aachen, Germany 

i^University College London, London WCIE 6BT, UK 

^^Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 
^^Department of Physics, University of Maryland, College Park, MD 20742, USA 
^^Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 
Canada 

^^University of Oregon, Department of Physics, Eugene OR 97403, USA 

20CLRC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 

Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 
^^Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
^^International Centre for Elementary Particle Physics and Department of Physics, Univer- 
sity of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 
^^Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
^^Universitat Hamburg/DESY, Institut fiir Experimentalphysik, Notkestrasse 85, D-22607 
Hamburg, Germany 

^^University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, 
Canada 

^^University of British Columbia, Department of Physics, Vancouver BC V6T IZl, Canada 

^^University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
^^Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hun- 
gary 

^"Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 
^^Ludwig-Maximilians-Universitat Miinchen, Sektion Physik, Am Coulombwall 1, D-85748 

Garching, Germany 

Max- Planck- Institute fiir Physik, Fohringcr Ring 6, D-80805 Miinchen, Germany 
^^Yale University, Department of Physics, New Haven, CT 06520, USA 



" and at TRIUMF, Vancouver, Canada V6T 2A3 
and Royal Society University Research Fellow 
and Institute of Nuclear Research, Debrecen, Hungary 
and Heisenberg Fellow 

^ and Department of Experimental Physics, Lajos Kossuth University, Debrecen, Hungary 
^ and MPI Miinchen 

^ and Research Institute for Particle and Nuclear Physics, Budapest, Hungary 



^ now at University of Liverpool, Dept of Physics, Liverpool L69 3BX, UK 

' and CERN, EP Div, 1211 Geneva 23 

^ and Universitaire Instelling Antwerpen, Physics Department, B-2610 Antwerpen, Belgium 
now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, USA 

' now at University of Toronto, Dept of Physics, Toronto, Canada 
current address Bergische Universitat, Wuppertal, Germany 

" and University of Mining and Metallurgy, Gracow, Poland 



1 Introduction 



The Standard Model (SM) has only one complex doublet of Higgs fields, resulting in one 
physical mass eigenstate, the neutral Higgs scalar boson [1]. Theoretical limitations of the 
SM have prompted the development of many other Higgs models. Possible extensions of 
the SM include the Two Higgs Doublet Models (2HDM). These models predict two complex 
doublets of scalar fields resulting in five physical Higgs bosons: two neutral CP-even scalars, 
h° and H*^ (with mh < mn), one CP-odd scalar, A°, and two charged scalars, H='= [2]. The 
Higgs sector of the Minimal Supcrsymmetric extension of the Standard Model (MSSM) is of 
2HDM(11) type where the introduction of supersymmetry adds new particles and constrains 
the parameter space of the model. Due to the Higgs boson self-coupling, decays of Higgs 
bosons into other Higgs bosons become possible if kinematically allowed. In the 2HDM(n), 
the couphng of the h° to A°A° is proportional to: 

[g mz/2cos(^w)]cos(2/5)sin(/3 + a), (1) 

where g is the standard SU(2) gauge coupling, mz the mass of the Z°, Oy^ is the weak mixing 
angle, a is the mixing angle that relates the CP-even Higgs states H° and h° to the field 
doublets, and tan /3 is the ratio of the vacuum expectation values of the Higgs scalar fields. 
When kinematically allowed, the decay h° — > A°A° may dominate. 

The Standard Model predicts that Higgs bosons accessible at LEP centre-of-mass energies 
decay preferentially into the heaviest available fermions since the coupling to the Higgs boson 
is proportional to the fermion mass. This fact has motivated the vast majority of the SM 
Higgs boson analyses to focus on Higgs boson decays via b quark and r lepton pairs [3]. 
OPAL has also performed flavour-independent searches to explore other possibilities [4]. 
Despite this effort, the region with uia < 10 GeV/c^ and between 45 and 86 GeV/c^ 
remains uncovered within the 2HDM(II) and MSSM parameter scans. The analysis described 
in this paper is dedicated to a narrow Higgs boson mass region, allowing for the selection 
of events with very specific kinematics. This guarantees a higher signal detection efficiency 
while achieving an excellent background rejection. This region has also been investigated by 
flavour-independent [4] and decay-mode independent [5] analyses but with lower sensitivity. 

The properties of the Higgs bosons in the MSSM can be studied in the framework of a 
constrained model with seven parameters: Msusyj M2, ^, Aq, tan/5, mA and rrig. Msusy is 
the sfermion mass and M2 is the SU(2) gaugino mass parameter, both at the electroweak 
scale. The parameter /i is the supcrsymmetric Higgs boson mass parameter, is the trilinear 
Higgs boson-squark coupling parameter, assumed to be the same for up-type squarks and 
down-type squarks, and rrig is the gluino mass. As an example, our results are interpreted in 
the MSSM no-mixing benchmark scenario [6], which assumes that there is no mixing between 
the scalar partners of the left-handed and right-handed top quarks, and is determined by 
the following choice of parameters: Msusy = 1 TeV/c^, M2 = 200 GeV/c^, = -200 GeV, 
the stop mixing parameter Xt = Aq — //cot/3 = 0, 0.4 < tan/3 < 50, 4 GeV/c^ < mA< 
1 TeV/c^ and rrig = 800 GeV/c^. In this scenario, the region for uia < 10 GeV/c^ and 
TBh within 45-86 GeV/c^ is not excluded by the OPAL data, and a smaller mass range 
with tua < 10 GeV/c^ and rrih around 70-86 GeV/c^ also remains unexcluded in the LEP 
combined data [7]. In this region, the A° boson is too light to decay into bb and would hence 
have escaped detection by the analyses using b-tagging. Furthermore, the existing flavour- 
independent analyses lack the necessary sensitivity to detect or exclude such possibilities. 
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The region for below 45 GeV/c^ is excluded by LEP 1 [4]. Although this search was 
originally motivated by the MSSM and 2HDM studies, our results can be extended to any 
model that predicts the existence of at least one scalar and one pseudoscalar Higgs boson in 
the mass range of interest and also within models where the physical Higgs bosons are not 
CP eigenstates, hke the CP violating MSSM [8]. For illustrative purposes, we will use the 
MSSM no-mixing benchmark scenario as a reference in the rest of this paper. 

The paper is organized as follows: Section 2 describes the data samples and the OPAL 
detector. Section 3 gives details of the Monte Carlo simulations and Section 4 describes 
the event selection. Section 5 covers systematic uncertainties and the results are given in 
Section 6. 



2 Data samples and the OPAL detector 

The data used for this analysis were collected during 1998-2000 at LEP in e"'"e~ collisions at 
centre-of-mass energies (Ecu) between 188 and 209 GeV. The data sample is divided into 
four subsamples, namely 188 < ^cm< 193 GeV (201.7 pb"^), 193 < Ecm< 198 (75.1 pb"^), 
198 < EcM< 203.5 GeV (115.2 pb"^) and 203.5 < £;cm< 209 GeV (206.7 pb"^). The analysis 
is performed separately on each sub-sample and the results are combined. The choice of the 
subgrouping of centre-of-mass energies is justified by the fact that the production cross- 
section for the signal does not change appreciably within each of the four subsamples since, 
for this analysis, we study a region far below the kinematic limit. 

The OPAL detector [9] has nearly complete solid angle coverage and hermeticity. The 
innermost detector of the central tracking is a high-resolution silicon microstrip vertex de- 
tector [10] which lies immediately outside of the beam pipe. Its coverage in polar angle^ 
is I cos^l < 0.93. The silicon microvertex detector is surrounded by a high precision vertex 
drift chamber, a large volume jet chamber, and 2;-chambers to measure the z coordinates of 
tracks, all in a uniform 0.435 T axial magnetic field. The lead-glass electromagnetic calorime- 
ter and the presampler are located outside the magnet coil. In combination with the forward 
detectors, namely the forward calorimeters, a forward ring of lead-scintillator modules (the 
"gamma catcher") [9], a forward scintillating tile counter (the "MIP plug") [11], and the 
silicon-tungsten luminometer [12], the calorimeters provide a geometrical acceptance down 
to 24 mrad from the beam direction. The silicon-tungsten luminometer measures the in- 
tegrated luminosity using Bhabha scattering at small angles [13]. The magnet return yoke 
is instrumented with streamer tubes and thin gap chambers for hadron calorimetry and is 
surrounded by several layers of muon chambers. 

Events are reconstructed from charged particle tracks and energy deposits ("clusters") in 
the electromagnetic and hadron calorimeters. The tracks and clusters must pass a set of 
quality requirements similar to those used in previous OPAL Higgs boson searches [14]. In 
calculating the total visible energies and momenta, Ey^^ and Pvis, of events and individual 
jets [15], corrections are applied to prevent the double counting of energy of tracks with 

^OPAL uses a right-handed coordinate system where the +z direction is along the electron beam and 

where +x points to the centre of the LEP ring. The polar angle, 9, is defined with respect to the +z direction 
and the azimuthal angle, (j), with respect to the horizontal, +x direction. 



associated clusters [16]. 



3 Monte Carlo simulation 

Monte Carlo samples for signal and background were generated at four different centrc-of- 
mass energies, namely 189, 196, 200 and 206 GcV, chosen to be close to the mean centre-of- 
mass energy in each data subsample described in Section 2. 

We study only h°Z° production since, in the parameter space region of interest for our 
analysis, its cross-section is about ten times larger than that for h^A" production in the 
MSSM. The h° is forced to decay into two A'' bosons, h° A^A*^, and each A*^ can decay into 
any of the following channels: cc, r+r^ and gg. Resonances arc not included in the simulation 
of A° decays. For example, in the MSSM no-mixing scenario, for 3.3 GeV/c^ < ttia < 9.5 
GeV/c^, the A° branching fractions into cc and r+r~ are 0.5-0.9 and 0.4-0.05, depending on 
the value of tan/9. Below the r+r~ threshold, the A° decays nearly exclusively into a gluon 
pair. Two different Z° decay modes are investigated: Z° uu and Z° — > with £ =c or 
jj,. For each of the Z° decay modes, the six final states obtained by all possible combinations 
of the A'^ decays to gg, cc and t+t~ have been analysed. In the no-mixing MSSM scenario 
below the production threshold for bb, these final states account for between 75% and 100% 
of the total decays of the A° boson [17]. The corresponding Feynman diagram is given in 
Figure 1. 



Figure 1: The Feynman diagram for the processes considered in this analysis. 

Monte Carlo samples were generated with mA=2, 4, 6, 9 and 11 GeV/c^ and for mh = 
45, 60, 70, 80, 86 GeV/c^ at each of the four centre-of-mass energies considered. For each 
[mA, nih] combination and each Z° decay channel studied, we produced 3000 events for each 
of the six possible final states using the HZHA2 [18] generator and the full OPAL detector 
simulation [19]. 

The branching fraction BR(h° — > A°A°) is relatively constant for uia in the range of 1 to 
11 GeV/c^ for a given value of ruh. The e''"e~ — > h°Z° production cross-section does not 
depend strongly on rrih in the range 45 < nih < 86 GeV/c^ but increases with increasing 
tan P values. 

Monte Carlo simulations are also used to study the various Standard Model background 
processes. The 2-fermion events, e+e~ — > qq, are simulated with the KK2f generator using 
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CEEX [20] for the modelling of the initial state radiation and PYTHIA 6.125 [21] for the frag- 
mentation and hadronisation processes. Bhabha events are generated with the BHWIDE [22] 
and TEEGG [23] generators, e+e" fi^fi^ and e+e~ — * r+r^ events are simulated with the 
KK2f generator using CEEX and ISR-FSR interference. The 4-fermion samples are gener- 
ated with grc4f [24] for the i£ii, iiqq, iuqq, ui/qq and qqqq channels, where i — e, /i, t. One 
2-photon sample generated at 200 GeV is used for centre-of-mass energies of 189, 196 and 
200 GeV and the resulting cross-sections normalised according to the centre-of-mass energy. 
An independent sample is used at 206 GeV. The various 2-photon processes are modelled 
by the Vermaseren [25], HERWIG [26], Phojet [27] and F2GEN [28] generators. Typically, 
at each centre-of-mass energy, the generated 4-fermion and 2-fermion samples are 30 times 
larger than the data sample, and the 2-photon sample is two to ten times larger than the 
data. 



4 Event selection 

In the low A" mass region covered by this search the separation between the decay products 
of the A° tends to be small, and they are generally reconstructed as a single jet. The final 
event topology consists of two jets recoiling against the Z° decay products. The invariant 
mass of a single jet reproduces the mass of the A° while the mass of the combined two-jet 
system reproduces the mass of the h°. A mass resolution between 0.5 and 3.0 GeV/c^ is 
achieved for niA, while for mh the resolution is about 8-16 GeV/c^, all depending on the 
Z° decay mode and the Higgs boson masses. 

We perform three separate analyses based on three different Z° decay channels: Z° — > uP, 
and e^e~. The event selection criteria are described in detail in the next two sections. 
After this initial selection, a likelihood variable is constructed for each channel to increase 
the sensitivity. This likelihood variable is used to perform a scan and set limits. 

4.1 The Z° vV channel 

Events in the process h'^Z'^ — > A°A°i/P are characterised by two jets recoiling against an 
invisible Z°. Each event is forced into a two-jet topology using the Durham algorithm [15] 
and a 1-constraint kinematic fit is performed, requiring energy and momentum conservation 
and forcing the missing mass to be equal to mi. Background from cosmic ray events is 
removed by requiring at least one track per event and applying the two cosmic-ray vetoes 
described in [29] and [30]. The visible invariant mass is also required to exceed 5 GeV/c^ to 
match a corresponding cut used in the generation for the 2-photon MC samples. This cut 
has no effect on this analysis since it is far from the region of interest for the signal. 

The selection criteria used to search for h^Z'^ — > P^P^vv are mainly based on event shape 
variables and reconstructed masses. The first four preselection cuts guarantee general data 
quality, confinement within the detector region and rejection of 2-photon background. The 
other more specific selection criteria aim mostly at rejecting the other backgrounds. The 
cuts used are described below: 



1. The fraction of total visible energy in the forward detectors must be less than 60%, 
the event transverse momentum measured w.r.t. the beam axis must exceed 3 GeV 
and the visible energy of the event must exceed 20% of the centre-of-mass energy 

2. Events are rejected if they have energy deposits exceeding 2 GeV, 5 GeV and 5 GeV 
in the forward calorimeter, the silicon-tungsten luminometer and the gamma catcher, 
respectively. This requirement rejects events with initial state radiation or particles 
escaping detection in the beam pipe. 

3. The fraction of good tracks as defined in [14] relative to the total number of tracks 
should be greater than 0.2. 

4. The polar angle, ^, of each jet must satisfy | cos^| < 0.9 to reject events with jets par- 
tially contained in the detector. This is also a very powerful cut to reduce background 
contributions from 2-photon and Z/7* events with large initial state radiation. 

5. The invariant mass of the di-jet system after the kinematic fit must be in the range 
30-95 GeV/c^. This mass should reproduce mh except in the case of A° t^t~ 
where there are missing neutrinos. 

6. The aplanarity^ must be in the range 0.0002-0.03 and the event shape variable C [31] 
must be less than 0.8. 

7. The polar angle of the missing momentum vector must satisfy | cos^| < 0.97. 

8. The invariant mass of the more energetic jet must be between 0.5 and 13.0 GeV/c^. 
(The mass of the reconstructed jet corresponds to rriA within the detector resolution. 
Since we perform a mass-independent search, we allow for a broad range of values for 
mx-) 

9. The invariant mass of the less energetic jet must be between 0.5 and 10.5 GeV/c^. 

The numbers of events passing cuts 1 to 9 are shown in Table 1 for all centre-of-mass energies 
combined. The numbers of events found in the data are compared to the various backgrounds 
after each cut. The qq and contributions to the 2-fermion background are listed in two 
separate columns. The numbers of events expected for a signal hypothesis of = 6 GeV/c^ 
and mh = 70 GeV/c^ in the MSSM no-mixing scenario and the corresponding efficiencies 
are also shown. 



4.1.1 Likelihood selection 

A discriminating variable is formed by combining information from the following four vari- 
ables into a likelihood based on the Projections and Correlations Method [32]: 

1. The event shape variable C . 

^The aplanarity is defined as IA3, where Ai are the eigenvalues [Ai > A2 > A3 with Ai + A2 + A3 = 1] of 
the sphericity tensor S°'^ = '}2iP'iP'^i I Yl,i |PiP> ^''^^ is related to the transverse momentum component out 
of the event plane. 



o 







1/1/ channel for £'cm= 189-209 GeV combined 


cut 


data 


total 
bgnd 


background sources 


signal 
(efficiency) 


qq 


e+i- 


4f 


27 


1 


21584 


20733.2 


9879.4 


4265.1 


2522.5 


4066.3 


73.7 (93.2%) 


2 


13241 


12776.5 


7793.0 


1688.2 


2044.5 


1250.8 


69.2 (87.6%) 


3 


12990 


12597.8 


7714.4 


1644.2 


2012.9 


1226.4 


68.8 (87.1%) 


4 


9015 


8876.9 


5281.1 


1300.5 


1617.0 


678.4 


68.8 (87.1%) 


5 


7837 


7809.4 


4966.2 


1003.1 


1332.5 


507.6 


59.4 (75.2%) 


6 


5227 


5338.8 


3991.0 


206.0 


776.2 


365.6 


50.9 (64.4%) 


7 


1571 


1547.5 


731.8 


142.3 


585.0 


88.4 


50.0 (63.3%) 


8 


854 


839.8 


429.1 


111.1 


237.8 


61.8 


48.6 (61.5%) 


9 


475 


470.5 


271.6 


57.7 


101.6 


39.6 


47.7 (60.4%) 


C > 0.88 


18 


14.9 


0.1 


1.0 


13.8 


0.0 


38.0 (48.1%) 



Table 1: The numbers of observed events together with the Monte Carlo expectation 
for the various sources of background given for the combined sample (189 < -Ecm ^ 
209 GeV) for the T? — > vv channel. The corresponding efficiencies are given within 
parentheses for one signal hypothesis (ttia — 6 GeV/c^, mh = 70 GeV/c^) in 
the MSSM no-mixing scenario. The various cuts are described in the text. The 
likelihood cut is only given for illustrative purposes. It is not used to set the limits 
in Section 6. 



2. Tfie acoplanarity^ angle of tfie two jets. 

3. The invariant mass of the more energetic jet. 

4. The invariant mass of the less energetic jet. 

These variables are shown for the data, the different sources of background and the refer- 
ence signal in Figure 2 after all preselection cuts given in Section 4.1 are apphed and for 
all centre-of-mass energies combined. To form the reference distributions for the signal, the 
distributions for sixteen [mA,^^!!] mass hypotheses obtained for mh= 60, 70, 80, 86 GeV/c^ 
and mA= 2, 4, 6, 9 GeV/c^ are summed after relative re-normalization according to the 
integrated luminosity of the data and their production cross-section. For the background, 
three reference distributions are used: the the 4-fermion, and the combined 2- photon 

and qq samples. Hence, the data are compared to four reference distributions: three for the 
Standard Model background and one single distribution for the signal. The likelihood func- 
tion used to derive our result is formed separately for each centre-of-mass energy considered. 
The distribution of the likelihood input variables for all centre-of-mass energies combined is 
shown in Figure 2 for illustrative purposes only. The efficiencies for each mass hypothesis 
are also determined separately at each of the [mA,77T.h] mass hypotheses and for each centre- 
of-mass energy. A small correction has to be applied to the efficiencies and backgrounds due 

•^Thc acoplanarity angle is the absolute value of 180° minus the opening angle between the two jets in 
the plane transverse to the beam direction. 
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to accelerator-related backgrounds in the forward detectors which are not simulated. From 
random beam crossing events the correction factors have been evaluated to be 3.1%, 3.6%, 
3.6% and 3.2% for ^ = 189, 196, 200 and 206 GeV, respectively. 

For illustrative purposes, a cut at 0.88 placed on the hkelihood variable would reject most 
background events and retain sufficient efficiency at all mass hypotheses considered. The 
same cut would be used for all data sets and is chosen by maximising the signal purity times 
the efficiency of our signal reference at all energies. After this cut, 18 events would be re- 
tained in data compared to 14.9 expected from SM backgrounds. The hkelihood distribution 
function is shown in Figure 5 (a) for the data, the Standard Model backgrounds and the 
reference signal for the 189-209 GeV data combined. The likelihood cut is not used to set 
limits in Section 6. The signal efficiency ranges from 38% to 75% for mp^ — 6 GeV/c^(see 
Figure 6(a)). It drops down between 21% and 53% for mx = H GeV/c^. 

4.2 The fi+jj,- and e+e" channels 

The leptonic Z° channels, namely h°Z° A°A°e+e~ and h°Z° A°A°/i+/i~ are also in- 
vestigated. These events are characterised by the presence of two jets with invariant mass 
compatible with the h" mass and a lepton pair with invariant mass close to the mass. 

The electron and muon analyses share the first three preselection cuts listed below. The 
selection starts with lepton identification. First the event is required to have one isolated 
lepton in association with two jets by applying the same selection criteria used in [33] for the 
isolation and identification of a lepton in qqii/ events from W"'"W~ decays. The selection is 
based on the probability for a track to be correctly associated with an isolated lepton. The 
probability is obtained with a likelihood method based on kinematic and lepton identification 
variables. No requirement on the number of tracks is made to avoid biasing the selection 
against low multiplicity events. Then the identification of two isolated leptons with the same 
flavour and opposite charge produced in association with two jets is required and the selected 
events are forced into a two jet configuration using the Durham algorithm without including 
the two best lepton candidates. 

The next two cuts arc appHcd to ensure confinement within the detector while the remaining 
selection criteria are optimised for background rejection and differ for the two analyses. This 
was necessary in order to reduce the large 2-photon background contribution in the electron 
channel. All cuts are optimised to maximise purity times efficiency for a mixture of all signal 
hypotheses. 

1. The isolation and identification of two oppositely charged leptons (e+e~ or /I'^ii'') in 
association with two jets. 

2. To avoid events having particles lost in the beampipe, both jets must have | cos^j < 
0.99. 

3. The visible energy must be greater than 0.78 of the centre-of-mass energy. This cut 
ensures the event is well-contained within the detector and rejects some of the 4-fermion 
background. 
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Muon channel 



4a. The invariant mass of the more energetic jet is required to be less than 25 GeV/c^. 
5a. The invariant mass of the less energetic jet has to be less than 15 GeV/c^. 

Electron channel 

4b. The invariant mass of the two leptons should be between 66 and 115 GeV/c^. 

5b. The number of charged tracks in each jet should be less than 10. This cut drastically 
reduces the 4-fermion background. 

6b. The invariant mass of the more energetic jet is required to be less than 36 GeV/c^. 

7b. The invariant mass of the less energetic jet must be less than 30 GeV/c^. 

8b. The angle between the two jets must exceed 1.6 rad to reduce the 2-photon background. 



muon channel for E'cm^ 


=189-209 GeV combined 


cut 


data 


total 


background sources 


signal 






bgnd 


2f 


4f 


27 


(efficiency) 


1 


56 


61.0 


0.7 


60.0 


0.1 


10.1 (77.4%) 


2 


55 


60.0 


0.6 


59.4 


0.0 


10.1 (77.4%) 


3 


44 


47.4 


0.2 


47.2 


0.0 


9.5 (72.8%) 


4a 


33 


37.6 


0.0 


37.6 


0.0 


9.1 (69.8%) 


5a 


27 


30.5 


0.0 


30.5 


0.0 


8.8 (67.5%) 


L > 0.56 


4 


3.6 


0.0 


3.6 


0.0 


7.8 (59.8%) 



Table 2: The numbers of observed events together with the Monte Carlo expectation 
for the various sources of background given for the combined sample (189 < Eqm < 
209 GeV) for the muon channel. The corresponding efficiencies are given within 
parentheses for one signal hypothesis (m\ = 6 GeV/c^, mh = 70 GeV/c^j in the 
MSSM no-mixing scenario. The various cuts are described in the text. The 2- 
fermion sample contains both qq and t^l~ events. The likelihood cut is only given 
for illustrative purposes and is not used to set the limits in Section 6. 

The number of events passing each of these cuts can be found in Tables 2 and 3 for the muon 
and electron channels, respectively. The numbers of events selected in the data are compared 
with the total background expected from the 4-fermion, 2-fcrmion and 2-photon samples after 
each cut. The number of events expected for a signal hypothesis of mh = 70 GeV/c^ and 
mA = 6 GeV/c^ in the MSSM no-mixing scenario is also shown. The 2-fermion background 
consists of qq and t'^t" events, but only qq events survive the preselection. While in the 
Z° — > ii'^iJ,~ channel the two-photon background is fully rejected by the preselection, in 
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electron channel for Eqm'- 


= 189-209 GeV combined 


cut 


data 


total 


background sources 


signal 






bgnd 


2f 


4f 


27 


(efficiency) 


1 


100 


103.6 


1.9 


80.4 


21.1 


9.9 (75.9%) 


2 


99 


99.5 


1.8 


79.3 


18.5 


9.7 (74.4%) 


3 


77 


79.5 


1.4 


62.0 


16.1 


9.2 (70.6%) 


4b 


35 


35.9 


0.4 


29.5 


6.1 


7.4 (56.7%) 


5b 


23 


19.7 


0.2 


15.2 


4.3 


7.1 (54.4%) 


6b 


21 


17.8 


0.1 


14.1 


3.5 


6.3 (48.3%) 


7b 


20 


16.2 


0.1 


13.1 


2.8 


6.1 (46.8%) 


8b 


19 


14.3 


0.1 


12.3 


1.8 


6.0(46.0%) 


C > 0.52 


4 


3.6 


0.0 


2.5 


1.1 


5.1 (39.1%) 



Table 3: The numbers of observed events together with the Monte Carlo expectation 
for the various sources of background given for the combined sample (189 < -Ecm ^ 
209 GeV) for the electron channel. The corresponding efficiencies are given within 
parentheses for one signal hypothesis (ttia — 6 GeV/c^, — 70 GeV/c^) in the 
MSSM no-mixing scenario. The various cuts are described in the text. The 2- 
fermion sample contains both qq and l^C~ events. The likelihood cut is only given 
for illustrative purposes and is not used to set the limits in Section 6. 



the Z° — > e''"e~ channel some hadronic tagged two-photon events survive the preselection 
cuts. To decrease the large statistical error on this background, all Monte Carlo generated 
2-photon events have been used at each centre-of-mass energy. 

4.2.1 Likelihood selection 

A discriminating variable is formed as for the Z° — > channel by combining information 
from the variables listed below into a likelihood based on the Projections and Correlations 
method [32]. The sixteen signal hypotheses obtained for mh= 60, 70, 80, 86 GeV/c^ and 
mA= 2, 4, 6, 9 GeV/c^ are combined to form one single reference signal distribution for 
each input variable as described in Section 4.1.1. This is done separately for the two leptonic 
channels. The variables used as inputs for the two likelihood functions are described here: 



Two-photon events in which one or both scattered electrons are detected [34] . 



Muon channel 



1. Angle between the more energetic muon and the nearest jet. 

2. Angle between the less energetic muon and the nearest jet. 

3. Reconstructed invariant mass of the more energetic jet. 

4. Reconstructed invariant mass of the less energetic jet. 

Electron channel 

1. Invariant mass of the electron pair. 

2. Angle between the less energetic electron and the nearest jet. 

3. Reconstructed invariant mass of the more energetic jet. 

4. Reconstructed invariant mass of the less energetic jet. 

5. Angle between the two jets. 

The distributions of the input variables for the data, the total background and the reference 
signal after applying all preselection cuts are shown in Figures 3 and 4 for the muon and 
electron channels, respectively. For the background, only one reference distribution is used 
in both channels consisting of the 4-fermion background sample: too few 2-photon and 
2-fermion events survive the preselection in the electron channel to use them as reference 
histograms. The likelihood distribution function is shown for the data, the Standard Model 
backgrounds and the reference signal in Figures 5 (b) and (c) for the muon and electron 
analyses, respectively. 

For illustrative purposes, a cut on the likelihood value could be set at 0.56 for the muon and 
0.52 for the electron analysis to optimise background rejection and signal detection efficiency 
at all mass hypotheses considered by maximising the purity times signal efficiency of our 
signal reference. This cut would retain four events in data compared to 3.6 expected from 
SM backgrounds both in the muon and electron channels. This cut is not used to set limits 
in Section 6. For the muon channel, the signal efficiency ranges from 32% to 77% for = 6 
GeV/c^ (see Figure 6(c)) and between 29% and 75% for mx — H GeV/c^. For the electron 
channel, the signal efficiency ranges between 14 and 57% at — 6 GeV/c^ (see Figure 
6(e)) and between 4 and 46% at mx — 11 GeV/c^. 



5 Systematic uncertainties 

Many effects related to possible inadequacies in the simulation of physical quantities in the 
Monte Carlo samples contribute to the systematic uncertainty. These contributions are 
hsted in Table 4 for the Z° — > i/V, Z° — > /x^/x" and Z° — > e''"e~ channels. They are added in 
quadrature to obtain the total systematic uncertainty for each channel. The evaluation of 
each source considered is described here: 

1 Q 







-> i/v channel 


Z"^ 


H'^li channel 


Z"^ 


e+e channel 


source 


bgnd 


signal 


bgnd 


signal 


bgnd 


signal 


tan A 


±0.1 


±(0.0 - 1.3) 


±0.3 


±(0.0-1.6) 


±0.4 


±(0.0-0.7) 


K 


±0.2 


±(0.1 - 1.1) 


±1.2 


±(0.2-7.0) 


±2.7 


±(0.3 - 2.2) 


inputs to £ 


±8.0 


±(0.7-8.3) 


±8.9 


±(1.0-4.2) 


±10.2 


±(2.0-8.3) 


SM cross-scc. 


±2.0 




±2.0 




±6.7 




MC statistics 


±4.8 


±(1.4-13.5) 


±10.0 


±(1.4 - 11.9) 


±11.7 


±(2.0 - 10.5) 


lepton ID 






±0.9 


±0.9 


±2.3 


±2.3 


total 


±9.8 


±(1.7-15.9) 


±13.6 


±(2.2 - 14.5) 


±17.3 


±(4.4 - 13.7) 



Tabic 4: Contributions to the systematic uncertainties for the vv, 1^^ 
and Z° e^e~ channels expressed in percent for the signal Monte Carlo and the 
Standard Model background, as described in the text. The error on the signal 
corresponds to the range of values obtained for the generated [mA,^hi values for 
MSSM Higgs bosons in the no-mixing scenario. 



• Simulation of likelihood input variables. Each likelihood input variable is re- 
scaled in the Monte Carlo so as to reproduce the mean and variance of the distribution 
seen in data. This scaling is done at the level of the preselection cuts described in 
Sections 4.1 and 4.2. After the illustrative likehhood cut, the difference in the num- 
ber of selected background events and signal efficiencies before and after re-scaling is 
used to derive the systematic uncertainty. The observed variations are then added in 
quadrature. 

• Detector tracking resolution. The tracking parameters for all Monte Carlo samples 

are smeared, varying the resolution by ±10% in tan A and ±6% in k separately to 
evaluate their contributions to the track reconstruction. Here, tan A = cot 9, where 6 
is the track polar angle and k, is the track curvature at the point of closest approach 
to the origin. The variation sizes are based on a comparison between data and Monte 
Carlo using muon pairs and Bhabha events. 

• Background cross-section determination. The overall 4-fermion cross-section un- 
certainty is assumed to be about ±2% [35] , reflecting differences in calculations of the 
W'^W" and Z^Z" cross-sections when comparing results from various generators. The 
same error is assigned to 2-fermion Standard Model background. For the electron chan- 
nel a 20% error on the cross-section for tagged hadronic two-photon events is assumed 
based on a recent OPAL measurement [36]. 

• Lepton identification. Lepton identification is performed as in [33]. Systematic 
errors are assigned to account for observed differences between data and Monte Carlo 
simulation in lepton identification and tracking efficiency. The lepton identification 
mismodelling is studied using "mixed events" . These events are formed by combining 
two kinds of events recorded at ^/s = 91 GeV: a Z*^ — > qq event is combined with 
half of a Z'^ ^ i^i^ event to simulate a W+W^ qq£u event. The systematic error 
is obtained by comparing data and Monte Carlo lepton identification efficiency in 



mixed events and is estimated to be 0.29% for electrons and 0.24% for muons. The 
second contribution to the uncertainty accounts for tracking losses. To determine this 
contribution, 1? e+e~ and T? fi^fi^ events in data recorded at = 91 GeV are 
compared to Monte Carlo events. The difference in the tracking efficiency is used to 
extract a 1.1% systematic error for electrons and a 0.4% error for muons. These two 
contributions are added in quadrature for each lepton, then doubled since our selection 
requires two such leptons per event, leading to a systematic uncertainty of 2.3% for 
electrons and 0.9% for muons. 

• Monte Carlo statistics. The numbers of events passing the preselection as well 
as the size of the Monte Carlo sample before preselection are used to determine the 
contribution from statistically limited samples based on binomial statistics. Table 4 
gives the contributions at ^/s = 189 GeV. In the limit calculation the systematic errors 
for each different final state and centre-of-mass energy are used. 

To illustrate the effect of the systematic uncertainties after the illustrative hkelihood cut 
described in Sections 4.1.1 and 4.2.1, a signal hypothesis of rriA = 6 GeV/c^ and rrih = 70 
GeV/c^ would contribute 38.0 ± 0.8 events in addition to the 14.9 ± 1.5 events expected 
from Standard Model backgrounds in the — > i/i/ channel after combining all Monte Carlo 
samples at 189, 196, 200 and 206 GeV. For the combined leptonic channels, the signal would 
contribute 12.9 ± 0.4 events in addition to the 7.2 ±1.1 expected events from backgrounds. 
There are 18 candidates selected in the data for the invisible channel and eight candidates 
in the leptonic channels as shown in Tables 1, 2 and 3. 



6 Results 

The analyses presented here are designed to explore the possibilities of a low mass Higgs 
boson, namely for rriA below the bb threshold. No significant excess of events is observed 
in either the invisible or leptonic Z° decay modes. Hence we set limits within two different 
scenarios: a model-independent scenario and the MSSM no-mixing benchmark parameter 
scenario. We obtain 95% confidence level (CL) exclusion limits using standard statistical pro- 
cedures based on the likelihood ratio technique [37] as applied in other OPAL publications. 
The likelihood variables described in Sections 4.1.1 and 4.2.1 are used as discriminating vari- 
ables for the limit calculation without cutting on this variable. For the signal, the hkelihood 
variable and the efficiency are calculated for each centre-of-mass energy, each [mA,?77.h] hy- 
pothesis and each final state. The search efficiency is computed from Monte Carlo samples 
produced for each A° pair decay channel (namely A*^A*^ — cccc, r^r^r^r^ , gggg, ccr+r", 
ccgg, T"'"T~gg), each [niA, rrih] hypothesis and each centre-of-mass energy considered in this 
study. For the rriA and rrih points located between the mass points where Monte Carlo 
samples were generated, the efficiencies, the shape of the likelihood distribution and the 
systematic errors are interpolated using a weighted mean of the relevant quantity for the 
four nearest [m,\- rrii^] mass points. The likelihood variable for backgrounds and data are 
determined separately for each centre-of-mass energy. The efficiencies are shown in Figure 
6 at £^cM = 189 GeV for each A°A° decay channel versus rriA and mh, both for the missing 
energy and the leptonic analyses. Similar behaviour is observed for any other choice of the 
masses and centre-of-mass energy. The efficiencies are calculated neglecting cases where the 
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A° would decay to resonances. The search is still sensitive to the A° decays to resonances 
since the resonant states decay preferentially into gg, t+t~ and cc [38]. 



6.1 Model-independent limits 

We calculate limits on the cross-section for the process e+e~ — > h'^Z". The limits can be 
extracted in terms of a scale factor that relates the cross-section for the production of 
h^Z'^, in any specific theoretical interpretation of our experimental search, to the Standard 
Model cross-sections: 

(ThOzO = S^(Th'4jZ0- (2) 

The h'^ — > A'^A'^ branching ratio is assumed to be 100%. The limits are extracted for 100% 
branching ratio of A^A*^ into cccc, gggg, t~^t~t^t~ , ccgg, ggr+r~ and ccr'^r^. For each 
of the six final states studied, Figure 7 shows the iso-contours of 95% CL exclusion for 
in the uia and mh mass plane with 2 < uia < H GeV/c^ and 45 GeV/c^< mh < 
86 GeV/c^. The scan is performed in 1 GeV/c^ steps in mh and in 0.5 GeV/c^ steps in ttt-a. 
The T^T~T~^T~ channel has the largest exclusion power despite the fact that the selection 
efficiency is slightly lower than in the other decay channels since the signal is better separated 
from the background. 



6.2 MSSM no-mixing scenario interpretation 

We scan the region with 2 < tua < H GeV/c^ and 45 GeV/c^< mh < 85 GeV/c^ in 
the rriA versus mh plane for the MSSM benchmark parameter scenario. The maximum 
theoretically allowed value for mh in this scenario is 85 GeV/c^ [6]. The scan procedure 
is the same as that of the OPAL MSSM parameter scan [39]. The expected number of 
events for the signal is adjusted so as to correspond to specific production cross-section and 
branching ratios for a particular point of the parameter space. The 95% CL expected and 
observed exclusion regions are shown in Figure 8. The region for 45 < mh < 82 GeV/c^ is 
excluded for 2 < ttia < 9.85 GcV/c^, i.e., up to the bb threshold where A'^ bb decays 
become dominant. For 82 < mh < 85 GeV/c^, the region is excluded for 2 < rriA < 9.5 
GeV/c^. The whole region below the bb threshold was expected to be excluded but is not 
due to the presence of candidates in the missing energy channel (see the third bin from the 
right in Figure 5 (a)). 



7 Conclusions 



We have searched for the process e+e" — > h°Z*^ with Z° decaying into uu, e+e~, ^'^fi~ and 
h° decaying into A'^A" with uia below the bb threshold. Six different decay modes for the 
A^A'^ system have been investigated: cccc, gggg, t'^t~t~^t~ , ccgg, ggr+r^and cct''"t~. No 
evidence for the presence of a signal has been found and exclusion limits have been derived 
both in a model-independent way and within the MSSM no-mixing benchmark scenario. 
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Large areas of the parameter space investigated have been excluded. In particular, in the 
MSSM no-mixing scenario, the whole rectangular area for 2 < ttt-a < 9.5 GeV/c^ and 45 < 
iTT-h < 85 GeV/c^ is excluded at 95% CL, while the expected exclusion area is 2 < < 9.8 
GeV/c^ and 45 < rrih < 85 GeV/c^. These limits are the best obtained so far in this region 
of parameter space, which has not previously been excluded. 
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Figure 2: The four input variables used for the hkehhood function of the Z" vv 
channel after all preselection cuts are applied. The contributions from the Standard 
Model backgrounds are added and normalised to the data integrated luminosity. 
The results are shown here for Ecm=189-209 GeV combined. The contribution 
from the reference signal is scaled up by a factor of five. 
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Figure 3: The four input variables used for the hkehhood function in the muon 
channel after all preselection cuts for Ecm= 189-209 GeV/c^ combined, where the 
labels 1 and 2 refer to the more and less energetic muon and jet, respectively. The 
only contribution from Standard Model backgrounds surviving the preselection, 
namely the 4-fermion sample, is compared to the data. The contribution from the 
reference signal is also shown. 





Figure 4: The Gve input variables used for the hkehhood function after all pres- 
election cuts in the electron channel for Ecm=189-209 GeV combined, where the 
labels 1 and 2 refer to the more and less energetic electron and jet, respectively. The 
contributions of Standard Model backgrounds surviving the preselection, namely 
the 4-fermion, 2-fermion and 2-photon samples, are compared to the data. The 
contribution from the reference signal is also shown. 
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Figure 5: The likelihood distribution functions in the (a) TP vv (b) —>■ 
and (c) Z° —>■ e"'"e~ channels are shown for the data, the Standard Model back- 
grounds and the reference signal (the mixture of all signal hypotheses) for 
EcM=189-209 GeV combined. The backgrounds are added and normalised to the 
data integrated luminosity. 
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Figure 6: Signal selection efRciencies at Eqm = 189 GeV versus 1711^(171 a) forrriA = 6 
GeV / = 70 GeV / c^) for the missing energy and the leptonic channels. The 
efRciencies are shown without any cut on the likelihood variable for all six decay 
channels of a A^A^ pair with A^ decaying into cc, r^r^ or gg. 
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Figure 7: Upper limits at 95% CL for in the ttj-a versus rriYi plane, assuming 
100% decays of h° into A°A° and 100% decays of A°A° into (a) cccc (b) gggg (c) 
T^T^T^T~ (d) T~^T~ gg (e) ccr+r~ and (f) ccgg. The iso-contour lines are for values 
of < 1, 0.8, 0.6, 0.5 , 0.4 and 0.2. These limits are derived using the combined 
results from T? —>■ vv, Z° SLnd TP —>■ e~^e~ channels and for centre- of- mass 

energies between 189 and 209 GeV. 
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Figure 8: Expected (dashed contour) and observed (light grey area) excluded re- 
gions at 95% CL in the ttia versus plane for the MSSM no-mixing benchmark 
scenario. These limits are derived using the combined results from T? — > vv, 
TP — ^ fi~^fi~ and T? — e^e~ channels and for centre- of- mass energies between 189 
and 209 GeV. The theoretically inaccessible regions and the region excluded by 
LEP 1 are also shown by darker 



